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Accelerated testing must address the failure mechanisms active within the devices undergoing tests in 
order to assess lifetimes in a meaningful way. The assumption of constant temperature, thermally 
activated lifetime, based upon the Arrhenius assumptions, does not always provide the necessary 
understanding to interpret accelerated tests in microelectronics. Plastic encapsulants, dielectric polymers, 
and underfill materials are subject to delamination and cracking with thermal cycling. Crack propagation 
during use environment exposure, drives the potential for failure of microelectronic devices and is 
therefore a necessary focal point in qualification and life testing. This paper reviews the available research 
in the application of fracture mechanics to this class of problems in microelectronics including relevant 
test data. In addition, useful acceleration factor models are derived for polymer crack propagation based 
on principles of linear elastic fracture mechanics. Further, a simple approach to estimating the minimum 
temperature cycling ranges, necessary to propagate a crack, is also presented. Finally, a methodology of 
applying acceleration factors to develop testing plans is shown, with an example in spaceflight for a 
cubesat in low Earth orbit. Overall, this is a paper that shows a useful and appropriate process for creating 
physics of failure based life testing for delamination and cracking failures in microelectronic polymers in 
a temperature cycling environment. 


Introduction 


Polymers play a crucial role in enabling modern microelectronic devices, from encapsulating simple diodes 
to the most complex digital processors. Polymers are a mainstay of modern microelectronic packaging [1, 
2]. These materials perform the crucial function of providing structural enclosures that protect the active 
silicon element from the environment, while enabling handling and interconnection of the devices for 
assembly into functional circuits. The mechanical and dielectric properties of polymer materials, along 
with the range of economic options in processing, are what make them so important to the 
microelectronic industry. As shown in Figure 1, an advanced device will encompass a number of layers of 
materials and interfaces within the structure. Adhesion and integrity of the polymer system is critical to 
the reliability of these devices [1]. Due to mismatches in properties and contact morphology these 
material interfaces are vulnerable to failure through delamination and crack growth under various loading 
patterns. The crack initiation and growth rate follows the theories of fracture mechanics which can be 
used to directly understand the durability of the system. These principles are subsequently explained 
along with useful applications to testing of microelectronic devices that depend on polymer packaging. 
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Figure 1: Simplified schematic of over molded flip chip ball grid array. The inset shows typically layered 
structure on which the reliability of these devices depend. 


Function and Properties of Microelectronic Polymers 


The major classes of polymers used in microelectronic packaging include epoxies and polyimides. For 
many types of advanced devices, polyimides are applied, as a film, over the hard layers of inorganic 
passivation, as one of the last steps in fabrication of the wafer. This application provides for defense of 
the active circuits on the surface of the die, against subsequent handling. This enables packaging and 
allows for the assembly of the silicon device into useful products, similar to the structure shown in Figure 
1. In addition, polyimide provides for a barrier against moisture during handling, packaging and service. 
Polyimide films also reduce stresses imparted to the surface of the circuits by the surrounding 
encapsulations that make up the package structure. Variations in polymer chemistry which make up this 
class of polymer (polyimides) will not be discussed in detail, but it is important to note specific polymer 
compositions may have specific responses to loading and environment, given different thermal and 
mechanical properties that are a result of the chemical variations in the polymeric chains [2, 3]. 


Polyimide chemistry can be made to be photoimageable. This combined with superior dielectric 
properties, and thermal stability, enable high yield processing that is key to cost effective advanced 
devices. Polyimide is preferred for advanced packaging applications such as the addition of solder bumps 
for flip chip. They are also used in many other types of devices including high voltage applications, PIN 
diodes and microelectromechanical systems (MEMS), owing to processability, mechanical and dielectric 
compatibility, and thermal stability. 


Epoxies are also highly versatile and diverse in terms of composition, properties and processing, enabling 
their use for structural functions and final packaging or encapsulation. Epoxy based adhesives are used 
for structural bonding, including die attach for conventional wirebond devices, as well as underfill for 
advanced flip chip devices. Figure 1 shows a typical high volume, high input-ouput (IO) pin count package 
schematic common in microelectronics. This is an example of an encapsulated device in which interface 
adhesion plays a critical role in the reliability of the device. There are many other devices types ranging 
from optocouplers to MEMS that depend on polymer integrity to ensure reliability. 


Die bond applications for epoxy adhesives are well understood, as this has been a mainstay of packaging 
for many years. However, development of flip chip has driven the need for improvement in understanding 
of reliability considerations for epoxies used in underfill applications. In the case of adhesives, pre-formed 


solids or liquid applications are used. For advanced flip chip devices, underfills are often applied in liquid 
applications which depend on capillary action to fill the gap between the device and chip carrier [2]. Filling 
the gap reduces the opportunity for moisture to reach the active surface of the die and provides for the 
structural integrity that relieves plastic strain in the solder joints. 


The properties of epoxy adhesives and liquids are modified with fillers to provide control of thermal 
expansion properties. A variety of particulate fillers can be used to control thermal expansion to more 
closely match other components in the package structure. This is the key to providing reduction and 
control of stress, given differences in thermal expansion between the die and the chip carrier. For many 
types of devices, printed circuit board laminates are preferred for cost considerations, but ceramics are 
also employed for specialty applications. 
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Figure 2: Solder bump electron micrograph showing the bump and surrounding epoxy underfill. The inset 
shows the epoxy matrix heavily loaded with particulate fillers to control thermal expansion of the matrix. 
(Courtesy of J.O. Suh, NASA JPL.) 


Figure 2 shows a solder bump section revealing the adjacent underfill that is heavily loaded with 
particulate fillers for control of thermal expansion of the epoxy matrix. Fillers may also be used to impart 
electrical conductivity or to enhance thermal conductivity. These fillers also affect the mechanical 
behavior including fracture as the crack tip interacts with the particulates. 


Epoxy molding compounds are most often used for final packaging of advanced devices but liquid 
applications are also used for a variety of applications for final encapsulation (glob-top). Transfer molding 
is often used to facilitate rapid processing and yield. The composition of molding compounds or 
encapsulants can be substantially different from adhesives, but fillers are also commonly used to control 
their properties. 


Failures in Polymer Packaging 


The system of materials in the realm of modern polymer encapsulated packages creates complex 
multilayer structures. The reliability and lifetime for these devices depends on the integrity of these 
multilayer polymer structures. However, failures and defects can and do occur in application and testing. 


Delamination or failure of a bonded interface is a primary failure mechanism of packaged structures. 
Therefore, adhesion of the layers in structure is critical to the reliability of many classes of semiconductor 
devices. Delamination precedes other failure mechanisms in microelectronic devices, which activate as a 
result of loss of the functions that polymers provide in the devices. Delamination will generally mean a 
loss of the dielectric or electrical insulating function, loss of structural integrity and stress control and a 
loss of the moisture barrier. Hence, for flip chip devices, delamination precedes corrosion of circuit 
elements and fatigue of the solder bumps [1, 4, 5]. An example of failure of underfill in an unencapsulated 
flip chip device, assembled to a ceramic substrate and employing an epoxy underfill, is shown in Figure 3. 
The micrograph prominently shows the fillet loss. However, the propagation of the crack along the epoxy 
to chip carrier interface in this device is of greater concern, as the propagation of the crack will increase 
the likelihood of solder bump failure with the loss of structural integrity afforded by the underfill. 


Cracking within the polymers, as well as along interfaces, is also a failure mode of concern. For example, 
cracking of polyimide passivation will cause a loss of the passivation protection it affords the die [6] 
leading to electrical leakage currents and corrosion. In the case of some types of high voltage devices, this 
will precede a catastrophic electrical failure. Passivation cracking in a polyimide coated die stack is shown 
in Figure 4. 


[ vie ] Lifted Fillet 


Solder 
Bump 


Underfill 


Figure 3: Unencapsulated flip chip structure with failed underfill. (Courtesy of Scott Popelar, Cobham 
Semiconductor Solutions (formerly Aeroflex)) 


Figure 4: Cracking in a polyimide passivation on a die stack. Cracks are decorated to enhance visualization. 
(Courtesy of Ahmed Amin and Henning Leidecker NASA GSFC.) 


Thermal stresses are predominant in driving failure in polymer encapsulated electronics. Thermal stresses 
arise from differences in the coefficients of thermal expansion (CTE) of the materials and the changes in 
temperature that occur from use environments and operation of the devices. A typical portion of an 
underfill is depicted in Figure 5. The difference in the thermal expansion between the epoxy underfill and 
the silicon die can easily exceed 10 to 15 ppm/°C depending on the amount of fillers in the epoxy and the 
specific epoxy composition, structure and processing. 
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Figure 5: Multilayer structure with typical coefficients of expansion. 


In simplest terms, the stress range developed in any particular material that is part of any stacked 
structure similar to Figure 5, are expressed in equation 1. They are proportional to the effective 


modulus of elasticity, E = El, Say where E is the Young’s modulus of the material and v is the 


Poisson ratio, as well as the differences in CTE, Aa, of the materials in the structure and the 
temperature changes, AT, that the structure experiences. 


AS = EAaAT (1) 


The mechanically stiffer elements drive the forces developed in the individual materials, as determined 
by the modulus of elasticity and thickness. When the system is not symmetrical, bending stresses will tend 
to develop, introducing additional complexity. Approximate closed form solutions to stresses in a 
multilayer structure are given by Suhir [35, 36]. The actual state of stresses in a multilayered structure 
may be complex and material properties can be dependent on temperature and time, as well. Detailed 
analysis of a microelectronic package are often best handled in carefully prepared finite element models 
(this topic will not be further developed and reader is referred to the references) [6]. 


Elements of Fracture Mechanics for Polymers 


Linear Elastic Fracture Mechanics (LEFM) has been widely applied to polymers to describe the growth of 
cracks in cyclic loading of laminates [7] and in studying materials that comprise microelectronics interfaces 
[8, 9, 10, 11 and 12]. Some basic elements of LEFM are explained below as necessary for understanding 
application to accelerated testing of microelectronic devices. 


In general, crack growth may occur in three specific fundamental modes. The opening mode or Mode | 
occurs under tensile loading in which the opposing faces of the crack are separated with application of 
the load, as reflected in Figure 6. In general, it takes less energy to propagate a crack in Mode | than in 
other modes. Mode II growth is also known as the edge sliding mode. This occurs in shear and is often 
associated with cracking in microelectronic type structures, under thermal cycling or vibration loads. The 
third mode, Mode Ill, also occurs in shear and is known as the tearing mode. 


Mode III ta 


Figure 6: Primary modes of crack growth. 


Crack growth in polymers used in microelectronics becomes a complex process of mixed modes of loading 
and crack growth. While shear loading apparently dominates, structures are rarely symmetrical, which 
leads to mixed modes of loading and crack growth, combining shear and tensile loading. Tensile loading 
develops as a result of some bending that will occur in non-symmetrical stacked structures of materials 
with different coefficients of thermal expansion (CTE) of the type shown in Figure 5 [8]. 


The basics of crack development in polymers differs from metals. Even in shear dominated loading, crack 
tip behavior in polymers is dominated by tensile behavior and tensile fracture properties; hence, Mode | 
material properties are most the important in describing failure [7, 13]. This is due to the formation of 
secondary microcracks at the forefront of the primary crack tip. In Mode II loading of polymers, these 
secondary cracks are near an angle of 45°to the fracture surfaces of the primary opening crack. This clearly 
shows the dominance of tensile behavior, as this orientation reflects the influence of the principal 
stresses. Ligaments form between the secondary cracks, in front of the primary crack, and subsequently 
rupture under continuous loading. The coalescence of rupturing ligaments forms new primary crack 
surfaces. 


At the interface between two different polymers, such as polyimide and epoxy, crack propagation occurs 
adjacent to or within an interfacial or interphase layer that forms as a result of bonding the two polymers 
[5]. Failure does not necessarily occur with a clean separation of two separate materials in a properly 
formed bond. Snodgrass [9] showed that cyclic growth of cracks in the epoxy-polyimide systems, for 
example, favored crack growth along the interface in the epoxy side of the bond, through formation of 
ligaments in the epoxy. 


These fundamentals of polymer mechanics are clearly not dependent upon sliding surfaces in Mode II in 
the classical sense of LEFM. It is the premise of this paper, therefore, that tensile properties, which are 
easily measured, can be used to provide for information to formulate accelerated test designs, when 
evaluating complex layers of polymers in microelectronic structures. 


Theoretical Foundations - Derivation of Fracture Parameters 


Prediction of interfacial failure in an accurate manner is heavily dependent on appropriately determining 
the stress intensity factors or strain energy release rate, whichever is more appropriate for the particular 
case. These interfacial properties are more dependent on bonding/fabrication process conditions unlike 
homogeneous materials. This makes specific tests necessary to determine the interfacial fracture 
parameters, which are then derived through load-displacement measurements and applied to theoretical 
models. 


Stresses attain their highest value at the corner edge of a multimaterial junction in electronic packages. 
High stresses and inherent bond weakness due to manufacturing defects can result in crack initiations at 
the corners. One such region of great concern is the bond between silicon die and flip chip underfill 
material (reinforced polymer). Sharp three dimensional corners completely embedded in polymer are 
locations of highest stress in the package under thermomchanical loading. 


Fracture theories have their roots in elasticity solutions where stress singularities take the form 
7 A 
Oj & Tp (2) 


where rp is the radial distance from the corner to the point at which stress is measured. The power index, 
A, independent of the loading conditions, is a complex number and can be expressed as follows [14]: 


A=ytin (3) 


where y and n are the real and imaginary components respectively. The consequence of A being a complex 
number is the stress state becomes oscillatory in nature. For stress state inside the nonlinear fracture 
process zone close to the crack tip can be solved by letting ai = +m and a2 = -m [15] in Fig. 7 which is 


asymptotic in nature and can be expressed in the form 
1 Por ie|~ as 
A ae {Re[Krp'*]6;;(, €) + Im[Krp‘*|&j; (8, e)+ K16;j;' (8, €) (4) 


where rp and @ are the local polar coordinates for a coordinate system in planes perpendicular to the crack 
front. 
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Figure 7: Local coordinates (rp, &) for asymptotic interface crack solutions 


The stress in eqn (4) comes from the complex exponent of r that can be expressed as 


r'€ = cos(elnrp) + isin(elnrp) (5) 


The stress intensity factor K has real and imaginary parts Kiand Ki (K = K; + K;;) that are similar to the 
conventional Mode | and Mode II stress intensity factors in homogenous materials. For interfacial 
problems, K; and K, are dependent on each other and only the Mode III behavior is independent. The 
imaginary power r® depends on the oscillatory index € as follows: 


e = — Inf] (6) 


B here is knowns as one of the two Dundur’s parameters. Dundur’s parameters are unitless quantities 
defining the elastic mismatch and can be defined by- 

= Ha (k2+1)—pH2(kK2+1) 

Ho (ky +1)+p1(k2+1) 

B= Hy (K2—1)—H2 (K2—-1) 

bo (Ky +1)+py (Kz +1) 


(7) 


(8) 


where u is the shear modulus, k= 3-4v for plane strain or k= (3-v)/(1+v) for plane stress (v is Poisson’s 
ratio). The subscripts 1 and 2 depict the two different materials on the two sides of the universe. 


So let’s now look at an example of a bimaterial interface of epoxy polymer and silicon interface which is a 
common failure site in electronic packages. Using their standard elastic properties of Silicon (E = 128.3 
GPa, v = 0.28) and Epoxy (E = 12.0 GPa, v = 0.33) and incorporating them in eqns. (7), (8) and (9) we can 
calculate the Dundur’s parameters as a = 0.8235, B = 0.2044 and a relatively large oscillatory index of € = 
-0.066. At that level of €, there is a high degree of coupling between Mode | and Mode II stress intensity 
factors. 


There are several ways to calculate the complex stress intensity factor K by solving the boundary value 
problem of the configuration within the realms of theory of elasticity. 


Figure 8: Bimaterial interface with an interface crack 


Let’s say a pre-existing crack of length 2a exists between the two semi-infinite material systems with 
remote tensile stress being o and remote shear stress being t in both materials. The transverse stresses 
however does not affect the stress intensity factor. Near the crack, both materials are in inhomogenous 
state. This boundary value problem has been solved by Rice and Sih [16]- 


K = (1 + 2ie)Vma(2a)*(o + it) (9) 


Now we can express the remote load (T) as a function of its magnitude (co) and direction (w) as follows- 
o+it=Tel” (10) 

So the stress intensity factor becomes 
K = (1 + 2ie)Vma(2a) Tel (11) 

The amplitude of the stress intensity factor is defined by 
|K| = KR (12) 

The real-valued quantity has the familiar dimension: 


|K| = [stress] flength]1/2 (13) 
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Indeed, |K|relates to the energy release rate G as [17]- 


2*°E,  £E,~* cosh?(te) 

Using these theoretical equations, a good understanding of crack behavior can be predicted. However, 
for very complex and nonlinear geometries a more accurate prediction can be made through the 
development of a finite element model that captures the irregularities in geometries as well viscoelastic 
or viscoplastic material properties. [18, 19 and 20] are some novel attempts on developing a simple 
elastic-plastic-creep finite element analysis (FEA) model including the various surface roughness features. 
The trick is to perform a two-tier modeling with the global model capturing the entire structure with all 
its far-field loading and then the local model capturing the crack tip behavior based on the near-field 
stresses. The exact material layers involved maybe different but similar modeling framework can also be 
developed for silicon-polymer interfaces. 


Of course there are experimental ways to compliment this as well. Sandwich specimens can be developed 
to measure the fracture energy of an interface between a thin film of material 1 and a substrate of 
material 2. One can sandwich the film between two substrates of material 1 and use a layer of adhesive, 
such as epoxy to glue the bare substrate to the one covered with the film. The substrates are much thicker 
than the film and the adhesive, so the whole specimen is easy to load. The interfacial fracture energy can 
be measured if the crack runs on the desired interface. Another efficient way of measuring the interfacial 
fracture energy would be to perform a 4 point bend test. 


A combination of any of these techniques can lead to a good measurement of the fracture energy. The 
following sections discusses how the crack growth rate and therefore the life estimation can be deduced 
from the stress intensity factor (or strain energy release rate) value. 


Cyclic Environmental Conditions and Crack Growth Rates 


As described in the previous section, the distribution of elastic stresses near the crack tip in an elastic 
body are described in terms of the stress intensity, K. With knowledge of this parameter, the elastic stress 
components can be estimated near the crack tip. However, plasticity does occur in a local region of the 
crack tip, as reflected in simplified diagram in Figure 9. The radius of the plastic region depends on the 
yield strength and may be approximated as shown in equation (15) [21]. 


Figure 9: Simplified view of plastic region near the tip of a propagating crack. vi 


Where 7, is the radius of the plastic zone and S, is the yield strength. 


1 [x]? 
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Values for K can be calculated for different loading and geometric configurations from consideration of 
the theory of elasticity. In general, K has the units MPa m’and takes the following form: 


K=wW-Syna (16) 


where S is the nominal stress developed in the component under load and a is the crack length for an 
edge crack. w is a function of the ratio of the crack length, a, to the cross sectional width, W. W is the 
length that the crack must traverse to cause failure. ‘“%Y depends on the geometry of the component. 
Deriving functional forms of K for complex, layered microelectronic package structures is a complex 
modeling task. 


Under cycling conditions that induce stable fatigue crack growth, the stress range may defined as AS 
where AS=Smax-Smin. Correspondingly, the stress intensity range becomes 


AK = a- ASVtta (17) 


The Paris Equation describes the rate of crack growth per loading cycle, = , over the range of AS in terms 


of the stress intensity range, AK. This takes the form of equation (18) 
 _ A. [AK]™ (18) 
dN 
where A and mare constants and m is unique to a given material or interface. 


It is often easier to estimate strain energy release, G, in testing basic materials or interface behavior. In 
this case the Paris equation takes the form: 
14 C“ (AGT? (19) 
dN 
G is a function of K? depending upon geometry and materials. 


Figure 10 shows the general behavior of crack growth in materials for components that undergo cyclic 
loading. Below a threshold stress intensity, AKin, there is no crack growth; the crack remains quiescent. 
The linear Paris region is a region in which the Paris equation applies. This is a predictable and stable crack 
growth often termed sub-critical crack growth. Beyond this region, at the Fracture Toughness, Kc, crack 
growth is rapid and unstable. This behavior will be exhibited for each loading condition and crack opening 
mode reflected in Figure 6. 


Given Mode | loading, with plane strain conditions, Ki. is known as the Mode | fracture toughness or 
simply fracture toughness; it is most often characterized in testing for materials used in structural 
applications. Applying equation (16) under the condition K= Ky, allows for calculation of critical crack 
sizes, a, in which impending failure will occur in components under load. Predictions of life or cycles to 
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failure require the integration of equation (18) between an initial defect size and a critical crack size at 
which failure occurs. 
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Figure 10: Behavior of crack growth in materials under cyclic loading conditions. The Paris Equation applies 
over the linear region between the threshold and fracture toughness and extend over several orders of 
magnitude in crack growth. 


Acceleration Factors for Delamination 


While stress and subsequently life prediction is complex, a useful acceleration factor can be derived for 
crack growth and applied to thermal cycling conditions to support accelerated tests for microelectronic 
devices; acceleration factors deal with proportional quantities and not predictions. The Acceleration 
Factor, AF, in terms of the cycles, N, accumulated under test and use conditions, can be expressed as 
follows: 

AF = Suse (20) 


Ntest 
Combining equations (17) and (18) and integrating for N: 


A 1 1 
N= ach a ae da (21) 


Since A, W, a and rare all constants, hence the following can be stated- 


1 
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From equation (1) the nominal stress range is proportional to the temperature range of the temperature 
cycle when a body is undergoing temperature changes: 


AS « AT (23) 


Then it follows from equation (20) that for two different conditions, use and test: 


AF = Nuse = [z==)" (24) 
Ntest ATuse 


When strain energy release rate is the primary consideration: 
AG « (AK)? (25) 


Then, given equations (1), (19) and (25), the AF is expressed as: 


— Nuse _ = [*teese]™ (26) 
Ntest ATuse 

To make these AF factors useful for thermal cycling tests of microelectronic devices, m or n must be 

known. These values may be extracted from controlled materials fracture tests with specialized specimen 

configurations that are conducted to characterize the material interfaces. A commonly used specimen is 

a double cantilever beam (DCB) from which the strain energy release rate and crack length during testing 

can readily be determined. It should be noted that these tests are normally isothermal. In addition, 


humidity can play a significant role in crack growth. 


As an example of an isothermal crack growth test, Zhu et. al [22] tested a bilayer structure of a typical 
microelectronic polyimide film composition over SiN. The plot of the data are shown in Figure 11. 
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Figure 11: Strain Energy versus crack growth rate for polyimide over SiN using a DCB test. The 
data reflect a slope of n=3.2 in the stable crack growth regime. With permission of IEEE. 
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Applying equation (26), the acceleration factor AF is expressed as: 


AF = Suse = [stesse]" (27) 
Ntest ATuse 

Snodgrass et. al. [9] evaluated epoxy underfill to polyimide interfaces. Applying equation (19) resulted in 

a value of n = 3.9. For the materials tested, crack growth rates are apparently higher at the same energy 

release rate for the epoxy-polyimide interface then for the polyimide-SiN passivation layer interface. The 

acceleration factor for thermal cycling, derived as discussed is: 


— Nuse _ [ATtest 18 
= Neest ie | 8) 


The large value of the exponents indicates that life tests on actual microelectronic package elements can 
be conducted in a relatively short time, over temperature ranges compatible with polymers, to assure 
that crack propagation will not present a risk in the field. This is significant since crack propagation 
precedes other failure mechanisms. This also indicates some caution is warranted in selecting test 
temperature ranges for plastic encapsulated microelectronic devices to ensure that unrealistic over 
acceleration does not occur. This, in turn, may reflect the need for procedures that differ from ceramic 
enclosures for conducting qualification tests. 


Table 1: Epoxy material on SiN under different loading conditions. (Courtesy of Reinhold Dauskardt, 
Stanford University) 


a Slope (m) 
0.5 6.8 
0.3 7.7 
0.1 8.1 


Table 1 shows tests done for an epoxy material on SiN by the Dauskardt Group in Stanford University. R 
denotes the stress ratio for the tests [21]. 


R= sa (29) 


Smax 
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A comparison of the slopes of the curves in the Paris region indicate similar behavior across the range of 
loading conditions with a relatively small trend. 


Table 2: Benzocyclobutene (BCB) on SiO2 over different humidity ranges. (Courtesty of Reinhold Dauskardt, 
Stanford University.) 


Estimated 
Humidit Threshold 
ievel ian Strain Energy Slope (m) 
Release Rate 
(J/m7) 

1% 18.2 0.69 
30% 23.2 0.34 
40% 26.7 0.22 
60% 34.3 0.21 
90% 45.3 0.24 


Table 2 shows tests conducted on Benzocyclobutene (BCB) -SiO2, also by Dauskardt Group, over a range 
of humidity levels. It is clear that increasing humidity levels significantly shifts the threshold for which 
crack propagation is necessary for this system of materials. Further inspection also showed as the 
humidity levels increase above 40%, the slope of Paris region of these tests reflect very little change. 


The strength and fracture properties of a bond depend upon the interphase layer thickness [5] and can 
be affected by aging and humidity depending on the hydrophilicity of the materials in the system [4]. The 
behaviors, reflected in Table 1 and 2, indicate that acceleration factors derived from DCB materials tests 
may have wide applicability over loading and environmental conditions due to relatively small changes in 
the slope m over a wide range of loading and environmental conditions. 


Acceleration Factor for Polyimide Cracking 


Notomi et. al. [23] also explored cracking of polyimide using a specialized type of small compact specimen 
and applying conventional crack propagation measurement methods. Compact specimen designs are 
often used for evaluating fracture behavior of bulk materials. Applying the Paris equation, they reported 
values for m for polyimide materials of varying molecular weight. They reported values of 3.3 < m< 5.2 
and indicated crack growth rates decreased with increasing molecular weight. 


Amagai [6] developed a model for life prediction of polyimide passivation coatings under specific 
conditions of an impinging plastic strain field on the polyimide coating due to the presence of a 
propagating delamination under temperature cycling. This model was developed for a specific type of 
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plastic package known as Small Outline J-Lead with Lead on Chip (LOC). The Amagai model was developed 
for a range of polyimide properties using simulation and Statistical Experimental Design. 


The Amagai model, can be used to predict the number of cycles to produce a failing crack in a polyimide 
film, N. It is is expressed in terms of the nominal stress range developed in the polyimide, AS, in 
temperature cycling, the estimated radius of the plastic zone, rp, at the delaminating crack tip, polyimide 
thickness, t, and polyimide film tensile toughness, ©. Amagai applied finite element models for stress 
prediction. 


The Amagai model is expressed as follows for the specific polyimide tested: 
(AS x iyo ‘ N9-4497 =t:Q (30) 


While the model is derived for specific conditions, it yields a more general acceleration factor for cracking 
of polyimide film. 


Given equation (20) and (30): 


0.4497 
N me AS1.095 


Nie Ag 02812 (31) 


The exponent derived for the acceleration factor based on the Amagai model for cracking of polyimide 
passivation in thermal cycling is within the range of values reported for isothermal fracture tests reported 
by Notomi et. al. [23] 


Threshold Limits 


The threshold temperature range at which crack propagation will occur is of considerable interest when 
designing accelerated tests to evaluate risks of defects contributing to failure. An estimate can be derived 
from simple considerations providing for a lower bound of concern for the threshold temperature range. 
As shown in Figure 8, crack propagation does not occur below a threshold stress intensity AKi,. Felck et. 
al [24] reported the lower bound of the Mode | threshold stress intensity for engineering polymers is 
approximately AK;,>.0.1 MPa m*. Fleck et. al [24] provided several plots of various materials to bound 
cyclic properties for design purposes; A Ki=0.1 becomes a useful basis for estimating a threshold 
temperature range. 


In evaluating a lower bound, a model for stress intensity is needed. For purposes of analysis of threshold 
limits herein, a block like structure under Mode II loading served as a closed form example for evaluating 
a lower bound of interest on the temperature range; the model for compact shear specimens was used 
for this purpose. The stress intensity for a compact shear specimen was explored in detail by Jones and 
Chisholm [25] and Hoyniak [26]. A symmetrical compact specimen with a driving stress, S, to deform the 
specimen and induce Mode II loading, has a closed form stress intensity shown in equations (32): 
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K = 1.1S/ma (32) 


Considering thermal stresses as a driver, where delamination may occur between epoxy and a passivated 
silicon surface under cyclic loading, equation (1) can be applied as a first order approximation. Combining 
equation (1) and (32) yields the following: 


AK = 1.1: EAaAT: ya (33) 


Considering a polymer in which the properties under mixed mode crack propagation are dominated by 
tensile behavior of the materials, and allowing AK = AK, , the threshold AT for given crack sizes, a, can 
be estimated from equation (33). 


As an example, polyimide passivated silicon die stack encapsulated in non-reinforced epoxy, with a 
potential for delamination between the polyimide and epoxy, is shown in Figure 12. 


Polyimide Coated 
Active Si Structure 


Figure 12: Polymide passivated die stack encapsulated in non-reinforced epoxy. Delamination of length, 
a, can be detected by radiographic inspection. (Radiography courtesy of Ahmed Amin, NASA GSFC.) 


Considering that the stiffer elements drive stresses, the following approximations are made: Aa = 
Qepoxy — Asi and E= Eepoxy ; the threshold AT can then be evaluated for a delamination of length a, 
given AKy, = 0.1 MPa m?/2. Threshold AT values are plotted against various crack sizes for the structure 
in Figure 13, under the assumptions stated. 


AT (°C) 


a/W 


Figure 13: Estimated threshold AT at which delamination will grow, at the lower bound of threshold stress 
intensity for polymer materials, versus normalized crack length, for a block like silicon structure passivated 
with polyimide and encapsulated in unreinforced epoxy. 


Figure 13 reflects for the simple structure that the threshold values of AT necessary to propagate a 
delamination at AK,, = 0.1 MPa m1/2 is relatively small and within the range of many typical thermal 
environments to which electronic products might be exposed. This lower bound estimate for the 
temperature range is also well below most test and screening environments. 


Test Plans for Risk Informed Systems Qualification (RISQ) 
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Reliable and cost effective implementation of advanced packaging schemes is needed for most systems 
applications. This requires adequate testing and qualification strategies that can assure acceptable risks. 
Testing and qualification must be completed within schedule constraints to provide the necessary added 
value during development. The aforementioned acceleration factors can be used to facilitate this. 


Risk Informed Systems Qualification (RISQ) can provide a cost effective approach to maximizing the 
information obtained from testing system packaging designs, and for qualification, as well as addressing 
potential risks from defects that might arise in a given assembly process. The approach for RISQ Analysis 
is outlined in Figure 14. 


Based on failure 
raarexeatclalsiaass 


System 
Environment 


Determine Evaluate 
Acceleration Risk 
Factor Preference 


Package 
Architecture 


Figure 14: Risk Informed System Qualification (RISQ) Analysis process for developing 
accelerated testing for packaging. 


System environments and selected architectures are inputs to the analysis. The environment encountered 
by the application will govern the appropriate selection of test conditions and models needed to develop 
acceleration factors [27, 28] appropriate to the failure mechanism and test design. The architectures 
selected are important, as materials systems and potential defect types are also important considerations 
in testing and test condition selection. 


The conditions in which an assembly will be used depend upon several design factors of the system, as 
well as the external environment. Under the hood automotive environment may be substantially different 
from the cabin environment and for spacecraft, relevant system factors may include the position in which 
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the assembly is located on an orbiting vehicle, the size of the spacecraft, or the potential use of thermal 
controls. 


Qualification testing often cannot be run in a reasonable time frame to demonstrate actual desired 
reliabilities of devices. However, tests may be run to a set of risk preferences for a given project. Risk 
preferences for a given system design are held by the decision maker. However, these may be influenced 
by standardized organizational risk preferences. An example of risk likelihood classification is shown in 
Table 3 for unmanned spacecraft. As comparison, in automotive applications, the preferred test may be 
to demonstrate a probability of failure, p, less than or equal to 0.1 at a confidence level, CL, equal to 
95%. 


Table 3: Risk Likelihood Classifications Used in GSFC-STD-0002 [29]. Prescribed probabilities of failure for 
classification of safety events and technical or performance risks are shown. 


Safety Technical 


Likelihood | (Estimated likelihood of (Estimated likelihood of not 
safety event occurrence) meeting performance 
requirements) 


Seve | pom) | Po 


}4 High = | (102< P< 104) (25% < P, 50%) 
(103< Py < 102) (15% <P, $25%) 
[1 Verytow ] Pest) | ormer sz) | 


Using standardized risk preferences for qualification test design quickly informs projects concerning the 
impact of using new and advanced packaging technologies. In Table 3, ps denotes probability of failure 
for safety event and pr denotes probability of failure for precipitating loss of performance. 


Test scenarios can be developed using the well-known Binomial-Parametric Test design recommended by 
Kececioglu [30]. For this type of design, the desired component reliability to be demonstrated may be 
estimated from the Weibull distribution. The probability of passing the test at a given confidence is 
estimated from the Binomial Distribution. Hence, for a given sample size, z, the following model may be 
used: 
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r 
z—-k 


Neve TNs /8 
1-CL= (2 eg) |e n) (34) 
k=0 
where jy = Nuse 7 
(-InR(N)pemo)® 


In the above model, CL is the confidence limit, z is the number of samples, r is the allowable 
number of failures, Nr is the calculated number of cycles to be executed in test, without 
acceleration, to meet the required risk preference. B is the Weibull shape factor and Nuse is the 
required field life. R(N)pemo is the reliability to be demonstrated in the test as determined from 
the risk preferences for the organization or project where: 


R(N)pemo = 1—p (35) 


where p is the risk preference or tolerable likelihood of occurrence of an undesirable event. 


The AF can be redefined in terms of our RISQ process model as: 


AF = Ne. = [Stes] (36) 


Ntest ATuse 


Using Equations (34), (35) and (36), the CL, z, and Nrest can be traded off against each other for a fixed set 
of targets, Nuse and R(N)pemo, given a value for m. The value of m represents the failure mechanism the 
test is designed against. Finally, unless the test population is very large r=0. 


It is very useful to present an example in applying the RISQ process. In this case, an example of a cubesat 
application for low cost spaceflight is shown. Cubesats will normally employ plastic encapsulated 
microcircuits given cubesat cost, power, mass and performance requirements. Considering a low Earth 
orbit (LEO), the expected temperature variation for an internal printed wiring assembly in a small cubesat 
structure may be as significant as — 20°C to +30°C or, AT=50°C, with each orbit [31]. However, this will 
ultimately depend on the thermal design. While many cubesats are designed for short life, the expectation 
is to expand their capability, mission and lifetime. A 3-year mission life in LEO will require 17,520 cycles. 


In order to apply the RISQ process model for developing qualification schemes for the plastic packages for 
the cubesat example, the Weibull slope for failure due to crack growth in the polymers is needed. Rusanen 
and Lenkkeri [32] determined that the Weibull shape parameter for fracture of epoxy polymer flip chip 
bumps was 1.4 < B < 2.5. This range values may be similar to the expected Weibull shape parameter for 
fracture and crack propagation in polymer packages where failure is defined to be a crack of sufficient size 
to result in activating another failure mechanism or precipitating the expected failure mode. 


Selection of test conditions should be judiciously made in view of the principles of fracture mechanics of 
polymers. Unless sophisticated analytical modeling techniques are employed to interpret the test results, 
the maximum test temperature should not exceed the glass transition temperature. The CTE and modulus 
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of the polymers will dramatically change at this temperature, as the polymer undergoes a phase change. 
For application in which moisture and humidity is an operational consideration, the maximum 
temperature change can be kept below 100°C. This will afford the ability to control humidity levels during 
the test. In addition, elevated humidity may favor crack propagation and lower threshold stress 
intensities, Kin. For the cubesat example, the operational and test conditions selected are shown in Table 
4. 


Table 4: Operational and Test Conditions for a Cubesat Application 


Operational and Test Parameters Range or Value 
Maximum Orbital Temp Range (°C) -40 to +80 
Device Max Operating Limits (°C) -40 to +100 
Payload Operating Range (°C) -10 to +40 
Operational Cycle (mins) 90 
Payload Design Life (years) 3 
Payload Design Life (cycles) 17520 
Life Test Temp Range (°C) -40 to +80 
Half Cycle Hold Time (mins) 15 
Humidity Control (%RH) None 
Delamination Acceleration Factor (m=7.8) 923 
Polyimide Cracking Acceleration Factor (m=4.7) 61 


Table 5: Test Options for Cubesat Microelectronic Qualification 
for Delamination and Cracking in Polymers 


Thermocycling -40°C to +80°C for 3 Year Mission Life 


Sample Size, z 5 10 15 
Allowed Failures, r 0 0 0 
Weibull Shape Parameter, B 1.4 1.4 1.4 
Confidence Level , CL 90% 90% 90% 


Demonstrated Reliability , R(N)pemo 0.98 0.98 0.98 


Life at Reliability Level, Neiea (cycles) 17520 17520 17520 
Required Failure Free Test Time, Ne (cycles) 163462 99631 74579 
Delamination E -Polyimide Interf N 
elamination Epoxy-Polyimide Interface, Ntest ia 108 Bi 
(cycles) 
Polyimide Passivation Cracking, Ntest (cycles) 2680 1633 1223 


Options for a RISQ based life test qualification plan for plastic microcircuits are reflected in Table 5 fora 
risk preference of p=2%, which is the maximum allowable accepted likelihood for a failure event, in which 
safety is not the prime consideration, to remain in a very low risk regime. 


Table 5 reflects the large acceleration factors for delamination of the epoxy-polyimide interfaces for the 
qualification test design. This does not infer that the plastics are not capable of meeting the application 
requirement. It reflects the sensitivity of a propagating crack or delamination to the cyclic load. It can be 
expected, given this sensitivity, that defects would rapidly manifest themselves as growing cracks during 
thermal cycling. However, from experience, some satisfactorily processed plastic encapsulated devices 
can survive and provide electrical performance after severe thermal shock from -55°C to +125°C [33]. 
However, such tests are not necessarily quantitative. Testing, if not appropriately applied, may limit the 
range of packaging technologies that can successfully meet requirements of the application. 


In applying acceleration factors, it is important not to extend the test conditions beyond the model’s limits 
to explain the materials behavior. In general, for polymers subject to cyclic loading, it is best to select test 
designs in which the applied stresses are such that Nrest >10* cycles and Nuse <10° cycles. Acceleration 
factors which project significantly outside this range may not be accurate. Fatigue behavior should not be 
expected to be linear outside these bounds [34]. 


Since polymer failure is precursor to electrical failure, the devices under test must be inspected following 
the test. Qualification or life test runs should include periodic electrical functional tests, radiography, 
scanning acoustic microscopy, or destructive analysis to evaluate the devices as test cycles are 
accumulated. 


Conclusions 


The reliability of plastic encapsulated microelectronics depends upon the integrity of the polymers, 
primarily epoxies and polyimides, which enclose and protect the devices. Effective quantitative 
accelerated life testing must include consideration of the failure of these polymer constituents and their 
respective interfaces. This is achieved through consideration of the failure mechanism, crack propagation 
and growth, that dominate these materials and precedes electrical failures and other mechanical fatigue 
processes. The principles of Linear Elastic Fracture Mechanics (LEFM) have been shown to effectively 
describe these processes in polymers. 


The application of fracture mechanics, including the Paris Equation, yields quantitative, yet simple, 
acceleration factors of familiar form that can be used to bound appropriate quantitative accelerated life 
tests for microelectronic parts. These tests can demonstrate effective risk reduction to real systems where 
the parts are employed. The parameters of these factors can be extracted from fracture tests that directly 
characterize the materials used in the encapsulation of microelectronics. Furthermore, application of 
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LEFM reveals the sensitivity of crack growth in polymers to stresses, which can be capitalized on to build 
appropriate and reasonable quantitative testing schemes. 
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